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Optical  phonons  and  electronic  properties  in  double  heterostructures 
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,  ABSTRACT 

i 

Optical  phonon  modes  in  a  semiconductor  double  heterostructure  (DHS)  are 
examined  within  the  continuum  model.  The  interface  modes  found  here  can  account 
for  the  novel  phenomena  observed  in  right-angle  Raman  scattering.  The  Hamiltonian 
for  the  electron  interaction  with  phonon  eigenmodes  is  derived  and  employed  to 
study  polaronic  effects  in  the  DHS.  The  1  s  -  2p  + 1 r ans i t i on  energy  of  a 

magnetopolaron  bound  to  a  hydrogenic  impurity  in  the  quantum  well  is  calculated, 
and  excellent  agreement  with  experiments  is  obtained. 


1 .  INTRODUCTION 

It  is  by  now  well  known  that  in  polar  crystals  such  as  III-V  compounds,  the 
electron-optical -phonon  interaction  plays  an  important  role  in  determining  the 
electronic  properties  at  higher  temperatures.  Since  these  compounds  are  the  usual 
materials  used  in  the  fabrication  of  microstructures  such  as  heterostructures  and 
superlattices,  there  has  been  growing  interest  in  the  study  of  electron-phonon 
interactions  in  these  structures.  Until  recently,  the  bulk  Frohlich  Hamiltonian 
has  been  employed  in  almost  all  calculations,  even  though  the  presence  of  surfaces 
and  interfaces  is  expected  to  change  the  nature  of  lattice  vibration  modes 

dramatically. 1 ’ 2 

Although  surface  optical  (SO)  phonon  modes  have  been  discussed  some  time 

ago,3’4  they  have  not  been  included  until  fairly  recently  in  the  treatment  of 

polarons  confined  in  a  slab6  or  near  the  interface  in  semi- inf inite  systems.6’7 

In  recent  papers,8’®  the  optical  phonon  modes  in  a  double  heterostructure  (DHS)  of 
polar  crystals  have  been  solved  within  the  continuum  model.  It  is  found  that  there 
are  only  two  types  of  phonon  modes:  the  interface  (IN)  and  the  confined  bulk 
modes.  The  latter  may  be  either  longitudinal  (LO)  or  transversal  (TO)  phonon 
modes.  Therefore,  the  Frohlich  Hamiltonian  involving  only  the  bulk  LO  phonon 
cannot  be  applied  to  quantum  wells  and  other  microstructures . 

We  report  in  this  paper  the  Hamiltonian  derived  for  the  e  1  ec  tr  on  -  optical  - 
phonon  interaction  in  a  DHS.  Both  the  interface  and  bulk  LO  phonons  are  included. 
The  interaction  Hamiltonian  is  then  employed  to  study  the  phonon  effects  on 
electronic  and  optical  properties  in  the  DHS.  The  self-energy  of  a  free  polaron  is 
calculated  for  arbitrary  well  width.  It  approaches  the  bulk  value  of  the  side 
material  in  the  limit  of  d  -  0,  and  that  of  the  central  layer  in  the  limit  of  d  -» 
«.  For  a  magnetopolaron  bound  to  a  hydrogenic  impurity  in  the  GaAs/AlGaAs  quantum 


\j  . 

/ 


well,  the  transition  energy  calculated  from  this  Hamiltonian  exhibits  the  strange 
pinning  phenomena  in  strong  magnetic  fields  as  observed  experimentally. 


2.  OPTICAL  PHONON  MODES 

The  detailed  solution  of  phonon  eigenmodes  can  be  found  in  Refs.  8  and  9. 
Here  we  just  outline  what  is  essential  for  our  discussion  of  polaronic  effects. 
The  geometry  of  our  DHS  is  shown  in  Fig.  1.  Material  v  (u  -  1,2)  is  characterized 
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Figure  1.  Geometry  of  the  double  heterostructure. 


bv  the  dielectric  constants  «  (w) ,  «  and  e  ,  plasma  frequency  u>  ,  number  of  ion 

v  ou  o°°  pv 

pairs  n  per  unit  cell,  natural  frequency  WQ|/ .  and  polarizability  a The 

polarization  vector  P  is  related  to  the  relative  displacement  u(r)  of  the  ion  pair 
and  to  the  macroscopic  field  E 


p(r)  -  u(?) 


where  we  have  defined  n 

v 


[  1  +  a  n  ( A  -  A  )  ]  ^  and 
1  i/  v  ov  v 


2  2 

A  -  4»ru)  /u> 

i/  pi/ 


(1) 


(2a,  b) 


The  dielectric  susceptibility  x^(u)  is  given  by 

4iry  1  -  (A  )rf^  -  4jt/3  (3) 

V  OV  V  V 

and  is  related  to  the  dielectric  function  by  *^(w)  ~  f  («)  *  1-  With  the  frequency 
parameters  (2),  we  find  from  (3) 


£l/(w)  -  £«l/(wL  ■  w2)/(“L  ■ w2) 

(4a) 

e  -  1  +  4irn  a  /(I  -  ~~  n  a  ) 

aoi/  V  V  J  V  1/ 

(4b) 

The  electrostatic  potential  <f>( r)  that  gives  rise  to  the  field  E  satisfies 


E(r)  -  -V*(?)  .  (5) 

On  the  assumption  that  the  dielectric  function  e^(w)  does  not  depend  upon  the 
spatial  coordinates  within  each  layer,  we  have  from  (5)  and  V»D  -  0 

«„(* -  0  .  (6) 

Since  the  DHS  does  not  possess  translational  invariance  in  the  z-direction,  the 
two-dimensional  (2D)  vectors  k  and  p  are  introduced  in  the  wave  vector  k  -  (/c,q) 
and  in  the  position  vector  r  -  (p, z) .  The  potential  can  then  be  written  as 

r)  —  J  die  4>(k  ,  z)elK’ p  (7) 

Substituting  (7)  ir.  (6),  we  find  that  the  Fourier  components  of  the  potential 
satisfy  the  equation 

«  (w)  -  *2]  <fi(K,z)  -  0  ,  (8) 

vdz  > 

which  may  be  solved  by  assuming  either 

e  (w)  -  0  ,  v  -  1,2  (9a) 

v 


or 


Mr  -  *2]  4>0 c,z)  -  0 


(9b) 


It 

Because 


2  2 

is  seen  from  (4a)  that  (9a)  is  satisfied  by  the  LO  frequencies  w  -  . 

*  “L2*  v*-krations  on  both  sides  of  the  interface  have  different 


time  dependences.  The  boundary  conditions  that  both  4>  and  the  electric 
displacement  must  be  continuous  across  the  interfaces  then  require  the 

confinement  of  the  LO  modes.  On  the  other  hand,  Eq .  (9b)  yields  the  interface 
modes  of  vibration.  The  polarization  eigenvectors  of  all  these  modes  as  well  as 
the  dispersion  relations  have  been  worked  out  in  Ref.  8.  It  is  also  shown  there 
that  the  novel  slab  modes  observed  in  right-angle  Raman  scattering10  can  be 
interpreted  as  the  interface  modes.  As  discussed  in  Ref.  8,  in  the  vicinity  of  the 
center  of  the  first  Brillouin  zone,  the  antisymmetric  interface  modes,  which  are 
predominantly  transversal,  oscillate  at  longitudinal  frequency  while  the  symmetric 
interface  modes  are  longitudinal  but  oscillate  at  transversal  frequency.  Using 
Eqs.  (8)  and  (9),  we  find  after  a  straightforward  calculation  electrostatic 
potentials  corresponding  to  these  modes.  Thus  for  the  confined  LO  modes,  the 
Fourier  components  of  the  potential  are 


<f>  ( k,z )  -  4wC  (k) 
m  m 


sin(m?rz/d) 

< 

cos(m»rz/d) 
with  the  normalization  constant 


for  even  m 
for  odd  m 

2-W 


(10a) 


Icm(*)|  -  V2/d  [*2  +  (mw/d)  j 

Corresponding  to  the  interface  modes,  the  potentials  are 

fC  V"1*1  |z|  >  i/2 

<f>  ,(**z)  -  •  J 

0^(1  +  e2^el^  c°sh(*cz)  ,  |z|  <  d/2 

for  the  symmetric  modes  and 

\z\  >  d/2 


(10b) 


(11a) 


aJ 


-C  sgn(z)  e 
d  4(*.z)  -  4  aj 


■K  Z 


Caj(l  +  £2/£i^sinh^KZ^ 


z  |  <  d/2 


(lib) 


for  the  antisymmetric  modes.  The  function  sgn(z)  -  1  for  z  >  0  and  -1  for  z  <  0. 
The  index  j  -  1,2  in  Eq  .  (11)  labels  the  two  branches  of  the  symmetric  or 
antisymmetric  interface  modes  which  have  the  same  normalization  constant  as 
discussed  in  Ref.  8.  The  normalization  constants  are  given  by 


i  „  /"*\  i  hJL  *d/2  . 

IV01  -7k  ’  ' 


2  2 
tan( ““ )  +  2 

pi  p2 


IF 


(12a) 


Ica.oo 


4*  icd/2  .  .  ,£d.  .  r0  ^  ^2 

-  ^  e  ytanh(2  )  /  (2*>—  +  — 

w  ,  w  o 

pi  p2 


tanh(^)  \h  .  (12b) 


3.  ELECTRON -PHONON  INTERACTIONS 

To  derive  the  Hamiltonian  operator,  we  start  with  the  free  polarization  or 
the  free  phonon.  In  terms  of  the  polarization  vector  n  and  its  canonical  conjugate 
-♦ 

7r,  we  can  write  the  Hamiltonian  in  the  form 


V  ’  2  ^  I 


—  [**(r)jr  (r)  +  w  **(r)ir  (r>  ] 
Z  a  a  a  a 


n  and  ?r  can  be  regarded  as  quantum  mechanical  operators  and  can  be  expressed  in 
terms  of  their  Fourier  components  as 


?(r)  -  ^  yjrf/8* [am(ic)  +  a^(-*)  ]*m(r) 


•£/$</ 8  *Aw 


sj  (*)  [asj(,C)  +  Vj 


+  y>(/8>rAw  Ak)  [a  (*)  +  a*  .(-*)]»  (r)]> 

aJ  aJ  aJ  ) 


(14a) 


*(r)  -  ^  (am(«)  -  a^(ic))jrm(r) 


}  {l7K“sj  <: 


)/8?rA  [a  .( k )  -  a  .(-*)]*  .(r) 
m  Sj  sj 


+  ivtyw  .  (k)/8*A  [a  (k)  -  aj.(-*)]  n  (r)l 
aJ  aJ  aJ  aJ  ) 


(14b) 


In  Equations  (14),  we  have  introduced  the  area  A  of  the  interface,  and  the  creation 

^  — *  — *  — * 

(annihilation)  operators  a  («)  (a  ( k ))  for  the  confined  modes  and  a  .  .(«) 

mm  sj  ,  aj 

(a  .  .(*))  for  the  symmetric  and  antisymmetric  interface  modes  of  phonons.  They 

®  j  ’ 

obey  the  commutation  relations 


(15a) 


[aQ(/c)  ,a^(/c* )  ]  -  S(k-k') 


faf(K),  aj(2')]  -  (a  W,a.(i')]  -  0  . 

a  p  a  p 


(15b) 


Using  these  commutators  as  well  as  the  orthonormality  relation  for  the 
eigenvectors,  it  is  straightforward  to  show  that 


H  -  H  +  Htkt 
ph  LO  IN 


LO 


[a^(ic)a  (k)  + 
u/  mm  l 


(16a) 

(16b) 


k  ,m 


HIN  '  ^  +  ^,4  (*>  faL  ('c>a,4  (*>  + 


sj 


aj  aJ  aj 


2]} 


*>J 


(16c) 


The  interaction  energy  for  an  electron  interacting  with  the  polarization  at  a 
—¥ 

point  r  is  represented  by  the  Hamiltonian 


H 


e-ph 


e<*(r) 


(17) 


Regarded  as  a  quantum  mechanical  operator,  the  scalar  field  can  be  expressed,  with 
the  help  of  (1),  (4)  and  (5),  as 


4>( r)  -  \  y>t/ [am(*)  +  a^( -k)  ]4>„(r) 


m 


/c ,  m 


'S  («) 


+  )  -ty^/BxAw  («)  [a  .(*)  +  at  <-*)]*  .(r) 

^.1  SJ  SJ  SJ  SJ 


+  7>(/87rAw  .(«)  (a  .(*)  +  aL  (**)  *  (r)} 

aj  aj  aj  aj  j 


(18) 


Combining  (10a),  (11)  and  (18),  we  find  after  some  algebra  the  interaction 

Hamiltonian 


H  ,  -  h  __  +  H  TM 
e-ph  e-LO  e-IN 


(19) 


H 


e-LO 


i**p  f  \  n  /  \  ,mnz.  .  a  ,  a|.  -* 

s  i  )  B  ( k )  cos(”)  fa  ( k )  +  a  (-« 

^  £  m  dm  m 


)] 


m-1 , 3 . . . 


(20) 


B(«)  s  in(“~)  [  a  ( 
m  dm 


m-2 , 4 . . . 


*)  +  al('*)]} 


H 


e  -  IN 


i**p-/c(  |  z ) -d/2)  f  -  „  - 

e  ^.(k)  U„  ,(*)  +  al  ;(-/c) 


sj  s.j 


s.J 


-  sgn(z)B  (*c)[a  (*)  +  a*  (-«)]l  for  |z|  >  d/2  (21a) 

aJ  3  *  J  3  »  J  J 

's,3* 

«.  j 


-  B_,(/0  TT^T^T  [a.  ,(*)  +  ..(-«)]]-  for  |z|  <  d/2 


aj  sinh(«d/2)  a,j 


(21b) 


The  coefficients  can  be  calculated  from  (10b)  and  (12).  The  results  are 


0  2  4»re  >(g>t 

'v*)!  iv;>i  r--v  , 

Lis  Ad(*  +  ~ — )  L  ®i/  oi/J 

d 


1  _  _l_ 

€  £ 


(22a) 


|B  .  .  (#6)  |  - 

sj.aj 


- sJl-T-o‘,t<‘  |c  aJ(I> I2 

8>rAtiJ  .  .  (k)  J  >  J 

sj  .  aj 


(22b) 


which  leads  to 


lV:>i2  'S' '  *VJ>/1 7it‘"h(f>  +72i 


(23a) 


,Baj  («)  I2  -  XT  -  ^aj(*)/[£lcoth(r)  +  e2: 


(23b) 


For  convenience,  we  have  defined  in  (23)  the  dielectric  functions  e  ( u )  by 

V  J 

1  1  1 


£  (w)  £  (w)  -  £  £  ((d)  -  £ 

1/  IS  OIS  IS  ®l/ 


(24) 


4.  POLARON  IN  A  DHS 


We  now  proceed  to  calculate  the  polaron  self-energy  by  perturbation  theor 
The  unperturbed  ground- state  wave  function  is  taken  to  be 


-  ,  ik  -p 

|0q>  -  jk  ,0;  >  "  yx  e  f<z>  !n^  > 

K  ,  0  K  ,  0 


(25) 


where  electron  states  are  characterized  by  a  plane  wave  k  in  the  xy-plane  and  the 

—¥ 

ground  state  0  in  the  z-direction.  The  electron  wave  vector  is  defined  as  k^  - 
(k  ,k  )  and  the  electron  ground  state  wave  function  in  a  quantum  well  is  given  by 


f(z)  - 


f  C  cos(kzd/2)  exp ( ik^ ( | z | -d/2) ]  ,  jz|  >  d/2 


(26) 


[  C  cos(k  z) 


z |  <  d/2 


The  wave  numbers  k  and  k^  are  related  to  the  electron  subband  energy  by 


kz  -  Jl mA/»2  .  k'z  -  72»e<Vo. Ep/H2 


(27) 


with  E.  determined  by  the  transcendental  equation 


E ^  -  Vo  cos2(f  /2me E^2)  ,  i  -  1 ,  2 ,  .  .  . 


(28) 


The  normalization  constant  C  is  given  by 


C  -  /2k'/(k'd  +  2) 
z  z 


(29) 


The  phonon  state  |n  >  is  specified  by  n  phonons  with  wave  vector  k  -  (*,q)  and 

k ,  0 

index  labeling  the  phonon  mode,  which  may  be  one  of  the  interface  or  confined  LO 
modes . 

The  unperturbed  energy  of  the  electron  in  the  ground  state  (25)  can  be 
written  as 


(0) 


*k  Y  1 

-  ST  +  Ei  +  2  (n+2} 

e  u  * ,  /9 


(30) 


In  the  weak- coupling  approximation,  the  corrections  due  to  the  electron-phonon 
interaction  can  be  calculated  by  standard  perturbation  theory  as 


«  -  e<°> 

k  k 
e  e 


-j 

*.0 


<k  ,1 


H  ,  k  -* 

K.0  e'Ph  ..IS 


•  <t  >|2  /  [< 

„  a  v 


(0)  (0) 


-  e  -  Yu> 

k  -K  K,0J 


-  AE^  +  AEj  +  AE^ 


(31) 


where  we  have  assumed  that  trans is tions  to  or  from  a  state  with  more  than  one 
phonon  are  negligible.  Because  we  ars  only  interested  in  the  electron  self-energy 
due  to  electron-phonon  interactions,  matrix  elements  corresponding  to  interband 
transitions  are  all  ignored.  With  the  Hamiltonian  operators  given  by  (19) -(21), 
and  the  wave  function  by  (25),  we  find,  after  integration  over  the  phonon  momentum, 
that 
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The  parameters  a  and  7  are  related  to  the  Frohlich- type  coupling  constants 
involving  the  particular  phonon  modes  indicated  by  the  subscripts.  Explicitly,  the 
a's  are  given  by 
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and  we  have  defined  the  function 
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K  ^  and  measure  the  polaron  size  due  to  the  LO  phonon  in  materials  1  and  2, 
respectively,  while  a  ^  (a^)  denotes  the  usual  LO - phonon- e lec tron  coupling 
constant  of  the  Frolich  type  in  material  1  (2).  They  are  given  by 
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Substituting  Eqs .  (32)  into  (31)  we  find 
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Combining  (36)  and  the  unperturbed  energy  (30),  we  obtain  immediately  the  ground- 
state  energy  of  the  interacting  electron-phonon  system  as 
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where 
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is  the  polaron  effective  mass  and 
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Equation  (39)  is  now  employed  to  compute  the  binding  energy  of  a  polaron  in 
the  GaAs  quantum  well  of  a  GaAs/AiAs  double  heterostructure.  The  results  are 
plotted  as  a  function  of  the  well  width  d  in  Fig.  2  in  which  contributions  from  the 
confined  LO  modes  and  interface  modes  are  shown  by  dotted  and  dashed  lines, 
respectively. 
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Figure  2.  Polaron  binding  energy  as  a  function  of  the 
well  width.  The  dotted  and  dashed  lines  represent  the 
contributions  from  confined  LO  and  interface  modes, 
respectively. 


It  is  observed  that  for  small  d,  our  results  are  qualitatively  different  from 
what  can  be  found  in  the  literature.  In  fact,  this  is  the  first  calculation  that 
gives  in  the  limit  of  d  -  0  the  correct  self-energy  ,  the  corresponding  value  in 
material  2.  This  is  because  we  have  included  the  contribution  from  confined  LO 
phonons  in  material  2  in  our  Hamiltonian.  The  interface  modes  play  an  important 
role  in  the  small-d  regime.  Their  influence  quickly  diminishes  as  d  increases.  It 
is  important  to  note  that  interface  phonons  are  responsible  for  the  pinning  effect 
observed  in  a  quantum  well  on  the  transition  energy  of  the  polaron  bound  to  a 
hydrogenic  impurity.11  This  is  discussed  in  the  following  section. 

5.  TRANSITION  ENERGY  OF  A  BOUND  MAGNETOPOLARON 

Recent  experiments11'12  reveal  that  the  resonant  effect  of  the  electron- 
phonon-  interaction  on  the  ls-2p  transition  energy  of  a  magnetopolaron  bound  to  a 
hydrogenic  impurity  in  GaAs/AlGaAs  multiple-quantum-well  structures  exhibits  the 
pinning  phenomenon  at  an  energy  below  the  TO-phonon  energy.  To  our  knowledge,  the 
origin  of  this  surprising  effect  has  not  been  understood.13  Here  we  attempt  to 
calculate  this  transition  energy  in  the  GaAs  well  of  a  GaAs/AlGaAs  L'HS . 

Consider  a  donor  impurity  at  the  center  of  the  well  as  in  the  case  of  the 
sample  used  in  the  experiments.  A  magnetic  field  is  applied  along  the  axis  of  the 
quantum  well.  The  total  Hamiltonian  can  be  written  as 


(40) 
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represents  the  energy  of  a  hydrogenic  impurity  confined  in  the  square  well  given  by 
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the  free  phonon  energy  is  given  by  (16)  and  the  interaction  Hg  ^  by  ( 19) - (21) . 

The  unperturbed  energy  levels  of  the  hydrogenic  impurity  described  by  (41a) 
are  calculated  variationally .  For  the  trial  wave  function,  we  take 
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solution  to  the  square-well  problem  and  is  given  by 
,  z  <  d/2 

(43) 

,  z  >  d/2 


where  k ^  and  k^_  are  defined  in  (27)  with  representing  the  energy  of  the  first 

subband.  The  function  G(p,z,<f> )  describes  the  internal  states  of  the  donor  and  is 
chosen  to  reflect  the  symmetry  properties  of  the  system, 
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where  £  is  the  variational  parameter  to  be  determined  by  minimizing  the  ground 
state  energy. 


The  interaction  ^  is  then  treated  as  a  perturbation,  and  the  energy 

levels  are  calculated  by  means  of  Wigner-Brillouin  perturbation  theory.  The  result 


is 
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where  E^(B)  stands  for  the  corresponding  unperturbed  energy.  The  calculation  is 

rather  complicated,  and  detailed  discussion  of  the  procedures  will  be  published 
elsewhere.15  Here  we  just  mention  that  the  Hilbert  space  has  been  truncated  to 
only  the  lowest  three  states  in  our  second-order  perturbation  term.  In  other 
words,  only  Is,  2p+  and  2p  are  included.  For  the  particular  samples  used  in  these 


experiments , 11  ’  12 


all  other  states  including  2p^  are  pushed  up  by  the  quantum  well 

wav  too  high  because  of  the  broken  symmetry  and  cannot  be  reached  with  presently- 
available  magnetic  fields.  Furthermore,  interactions  of  the  electron  with  lattice 
vibrations  outside  the  well  under  consideration  are  ignored,  since  they  are 
important  only  for  a  very  narrow  (d  <  10  A)  well,  while  d  is  larger  than  100  A  in 

the  experiments  considered  here. 


Magnetic  Field  (I) 


Figure  3.  ls-2p+  transition  energy  as  a  function  of  the 

magnetic  field  B.  Data  are  taken  at  SUNY-Buffalo  and  at 
the  Francis  Bitter  National  Magnet  Laboratory  in  a  series 
of  measurements  as  indicated  at  the  lower  right  corner. 
The  heavy  solid  line  represent  results  from  the  present 
theory  and  the  straight  line  is  the  result  without 
electron-phonon  interaction. 


The  ls-2p+  transition  energy  calculated  by  computer  iteration  is  plotted  as  a 

function  of  the  field  B  in  Fig.  3  along  with  the  newly- improved  set  of  data.12  It 
is  remarkable  that  the  theoretical  curve  breaks  into  three  branches  separated  by 
two  gaps,  in  good  agreement  with  experiments.  The  occurrence  of  gaps  reflects  the 
energy  level  repelling  caused  by  the  strong  resonance  interaction  when  the 
transition  energy  matches  the  phonon  energy.  Since  the  electron  does  not  couple  to 
TO-phonons,  this  pinning  phenomenon  must  be  due  to  the  electron  interaction  with 
interface  phonons  which  oscillate  at  frequencies  between  the  bulk  and  The 


interface -phonon  coupling  alone,  however,  is  not  adequate  for  a  complete 
understanding  of  this  phenomenon.  There  are  still  small  discrepancies  between 
experimental  data  and  the  theory,  especially  in  the  lower  branch  near  the  gap. 
Minor  corrections  such  as  the  nonparabolic  effect  on  the  effective  mass  can 
probably  account  for  the  discrepancies .  This  and  other  possible  corrections  will 
be  discussed  in  a  forthcoming  article.15 
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